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INTERFACE  ENERGIES 


Cheng-Chih  Pei 


Department  of  Physics,  University  of  California  San  Diego 
La  Jolla,  California  92093 


Abstract 


We  developed  a practical  scheme  for  the  calculation  of  interface  energies.  It 
combines  the  theory  of  generalized  Wannier  functions,  the  generalized  recursion 
method  for  calculating  local  densities  of  states  and  electronic  density,  and  the 
local  density  functional  theory.  As  a first  application  of  this  method  we 
are  calculating  the  stacking  fault  energy  of  nickel  using  tight  binding  type 
Wannier  functions  for  the  d-electrons  and  ignoring  the  effect  of  the  s-electrons. 
The  d band  degeneracy  is  fully  taken  into  account.  The  method  also  allows  one 
to  handle  charge  transfer  effects:  a stacking  fault  perturbation  potential  is 
included  and  calculated  quas i -sel f-cons i stent ly.  Comparison  is  made  between  the 
moment  scheme,  the  non-self-consistent  scheme  and  the  self-consistent  scheme. 

We  find  that  the  self-consistent  scheme  affects  very  little  the  stacking  fault 
energy  compared  with  the  non-self-consistent  scheme. 


* Supported  in  part  by  the  Office  of  Naval  Research  and  the  National  Science 
Foundation. 


Chapter  I 


Introduction 

1 2 • 
Hohenberg  and  Kohn  and  Kohn  and  Sham  (KS)  developed 

a procedure  for  determining  the  ground  state  energy  of  a system  of 
interacting  electrons  subject  to  an  external  potential.  This  procedure 
involves  the  solution  of  a single  particle  Schrodinger  equation  for  an 
electron  moving  in  an  effective  potential  that  is  determined  self- 
consistently.  Traditional  methods  of  solving  this  one-particle 
Schrodinger  equation  for  crystalline  solids  have  involved  obtaining 
the  Bloch  band  structure.  A perfect  infinite  and  periodic  lattice  is 
an  essential  ingredient  in  this  method.  If  one  is  to  deal  with  non- 
periodic solids,  such  as  clean  or  overlayered  solid  surfaces,  inter- 
faces, problems  associated  with  bulk  impurities,  stacking  faults, 
vacancies,  etc.  the  eigenstates  and  eigenvalues  of  this  Schrodinger 
equation  are  very  difficult  or  impossible  to  solve  for.  In  addition  to 
this  difficulty  of  calculating  the  eigenstates  and  eigenvalues,  the 
iterative  procedure  described  by  KS  is  in  general  quite  time-consum- 
ing  to  implement  as  has  been  pointed  out  several  times.  In  this 
part,  we  want  to  develop  a practical  scheme  for  the  calculation  of  the 
ground  state  energy  for  transition  metals  in  a non-periodic  system. 

Because  of  the  localized  nature  of  the  cL  orbitals,  it  seems 
reasonable  to  study  the  role  of  the  d-electrons  in  the  properties  of 
transition  metals  in  the  tight  binding  approximation.  A recursion 
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method  to  calculate  approximately  the  local  electron  density  of  states 

and  charge  density  within  the  tight  binding  model  developed  by 
4 5 

Haydock  et  al.  and  Rehr  and  Pei  in  a non-periodic  system  without 
the  necessity  of  calculation  eigenstates  of  the  one-particle  Schrodinger 


equation. 


We  now  discuss  our  procedure  for  handling  the  iterative 


difficulty.  The  KS  self-consistency  loop  goes  as  follows: 

Starting  from  a trial  electron  density  /]fi0(£) , one  constructs 
the  effective  potential  and  finds  a new  A new  effective  potential 

is  then  generated  and  the  entire  procedure  is  iterated  until  self- 
consistency  is  reached. 

Since  the  ground  state  energy  has  the  stationary  property 
with  respect  to  /YH£)  , the  error  of  the  ground  state  energy  is  of  the 
second  order  in  the  error  of  'THf)  • That  means  even  without  any 
iteration,  we  still  can  have  very  accurate  results  of  the  ground  state 
energy  provided  we  choose  a good  initial  trial  00  (£)  . Again,  be- 
cause of  the  localized  nature  of  the  cA.  orbitals,  the  state  of  the  elec- 
tron in  a crystal  is  very  similar  to  that  in  the  free  atom.  Therefore, 
a superposition  of  partially  overlapping  charge  densities  for  individual 
atoms  placed  on  the  appropriate  lattice  sites  will  be  a very  natural 
and  good  choice  for  the  initial  trial  . 

A description  of  our  method  for  calculating  the  ground  state 
energy  is  presented  in  Chapter  II.  In  Chapter  III,  we  present  a 
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— 


simple  illustration  of  this  method  in  which  the  stacking  fault  energy 

10  0 

is  calculated  for  a nickel  in  the  3d  4s  configuration.  In  Chapter  IV, 
we  study  the  stacking  fault  energy  for  a nickel  crystal  in  the 
3d^*  ^4s^‘  ^ configuration  with  a self-consistent  method.  Finally, 
Chapter  V contains  a discussion  and  a comparison  with  some  related 
methods. 
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Chapter  II 


Calculational  Procedure 


The  general  scheme  of  our  method  is  indicated  by  the  "flow 
chart"  of  Fig.  1.  A more  detailed  discussion  of  this  flow  chart  is 
given  in  the  following  paragraphs. 

The  ground  state  atomic  wave  functions  are  taken  from  the 
Herman  and  Skillman  atomic  orbitals.  ^ These  functions  are  fitted  by 
a set  of  Gaussians  in  order  to  facilitate  the  subsequent  Hamiltonian 
matrix  element  computations.  The  atomic  charge  densities  of  elec- 
trons, described  by  these  atomic  wave  functions,  are  chosen  to  be 
spherically  symmetric  and  the  total  electronic  charge  density  in  the 
crystal  is  constructed  by  adding  the  overlapping  charge  distributions 
produced  by  neutral  atoms  placed  on  the  appropriate  lattice  sites. 

From  this  electron  charge  density  the  effective  potential  is 
constructed.  An  important  saving  in  computation  time  is  obtained  by 
separating  the  effective  potential  into  a sum  of  overlapping 

spherically  symmetric  atomic  effective  potential  \/gw  l£)  plus  a 
corresponding  correction  term  ^ (f)(this  term  comes  from  the 
nonlinearity  of  the  exchange  and  correlation  potential): 


V11-? 


(2.  1) 


whe  re 


r 


jfr * j ji + f-„c ( <2-2> 
AVJtt(n=/*(ni0tt))-Zl  ttxtl^a-A))  <2-3> 

'n“(i)=|,  . ,2-41 

Here  jL  runs  over  all  the  lattice  sites,  and 

yttxc U(£>)  = dl  ( 'H  £x tl"))/ 1 <2-5) 


where  &«W  is  the  exchange  and  correlation  energy  per  electron  of 
a uniform  gas  of  density  00  . The  first  and  second  terms  on  the 

right  hand  side  of  Eq.  (2.  2)  are  the  Coulomb  potential  and  the  third 
term  is  the  exchange  and  correlation  contribution  to  the  chemical 
potential  of  a uniform  gas  of  density  'A  . Again,  for  facilitating  the 
subsequent  Hamiltonian  matrix  elements  computation,  the  atomic 
effective  potential  is  fitted  to  the  following  'Gaussian"  form: 


(IV 


c,  a -^r> 

III  +u1Cie 


(2.  6) 


a , the  nonlinear  part  of  the  exchange  and  correla- 

tion potential  is  small  and  a weak  function  of  f , which  in  the  first 
order  approximation,  can  be  treated  a6  a constant  ^ A . This 
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constant  can  be  obtained  by  taking  an  average  of  the  values  of  ^ V ) 
over  a unit  cell. 

7 

At  the  same  time,  the  generalized  Wannier  functions  (GWF) 

of  dl  orbitals  are  constructed  from  these  normalized  atomic  orbitals. 

These  GWF  v£)  are  localized  respectively  about  the  lattice  sites 

I *.8 

£ and  orthogonalized  by  the  Liawdin  algorithm  with  respect  to  all 
the  orbitals  (including  Is,  Zs,  Zp,  3s,  3p,  and  3d  orbitals)  which  are 
localized  at  the  nearest  neighbors  of  . 

If  we  choose  these  orthonormal  GWF  as  localized  basis  func- 
tions, we  can  calculate  the  matrix  elements  of  the  Hamiltonian  H 
for  various  relative  orientation  of  sites  (see  Appendix  A).  In  our 
computations  only  the  nearest  neighbor  interaction  have  been  included. 

In  this  case,  the  matrix  elements  of  H can  be  expressed  in  terms  of 

the  self  energies  and  the  nearest  neighbor  hopping  integrals. 

4,9,10  11 

Several  techniques  ’ ’ ’ have  been  developed  to  calcu- 

late the  local  electronic  density  of  states  from  the  matrix  elements 
of  a tight  binding  Hamiltonian  in  a non-periodic  system.  In  the  follow- 
ing we  shall  use  a relatively  new  and  powerful  generalized  recursion 

4 

method  which  was  developed  by  Haydock  et  al.  and  generalized  by 
5 

Rehr  and  Pei  to  find  the  local  electronic  density  of  states  '71. (E)  and 
the  electron  charge  density  ,YH£)  . 

Before  proceeding  with  the  calculations  of  and  '*1(0  . 

let  us  briefly  review  the  basic  structure  of  the  recursion  method. 
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In  this  method  the  local  density  of  states  at  a given  site,  say  at  site  i.  , 
is  calculated  from  the  diagonal  Green's  function  matrix  element: 


7\ .(E)  = 2,  OV  (E) 

JL  * x 


Jo( 


/Y|?(E')  = fix™  - X Ivn  fl  (Eti-S) 


S-*o 


u 


(2.7) 


(2.  8) 


where  o(  is  the  index  describing  the  type  of  the  corresponding  orbital 
at  site  1 (for  a 5 -fold  d band,  there  are  five  types  of  orbitals  at 


£>)  i 


each  lattice  site),  and  where  '-tg  j 

/v 

fraction  representation: 

<ot» ilce-HT1! ->(,  t> 

_ L 


s evaluated  using  a continued 


E-v 


e-«r  - 


1,1 


(2.  9) 


Here  | denotes  a GWF  of  type  0^  at  site  1 . For  the 

moment  we  restrict  our  discussion  to  a particular  oL  , for  simpli- 
fying the  notation  the  explicit  o(  dependence  will  be  suppressed. 
The  coefficients  (Xj  - and  b.  - in  Eq.  (2.9)  are  precisely  the 

1,3 

coefficients  of  a 3-term  recursion  relation,  which  defines  a new 
basis  U,il  in  which  the  Hamiltonian  has  a "tridiagonal"  repre- 
sentation (a  tridiagonal  matrix  is  one  whose  non-zero  elements 
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appear  only  on  the  main  diagonal  and  the  two  sub-diagonals,  upper 


i 


! 

right  and  lower  left): 

where 

(IjlHlMl  , Ifl.ilHIU+iir 

%j  U, ili.il  ; hi  U, ill, ilUji-HlM+'l.  (Z.1D 

We  shall  refer  to  the  states  U.il  in  this  semi-infinite,  one- 
dimensional basis  as  "shell  states,  " ^ being  a shell  index,  since, 

roughly  speaking,  they  are  spatially  localized  about  the  shell  of 

neighbors  surrounding  the  atom  at  site  ^ . 


Following  the  generalized  recursion  method,  we  can  write 
the  single  particle  energy  6 and  the  electron  charge  density  01(X)  : 


M being  the  total  number  of  electrons  in  the  system  and 


^=2,  (Ur),4  *(r)  f p te) 'fl  te^e 


(2.  14) 


where  A,  ill),  the  normalized  shell  states,  can  be  expressed 

-•4 


as 


linear  combination  of  GWF : 


A ( r)  — ^ ' A>1 1 

hr  fiJW 


(2.  15) 


and  where  are  the  orthonormalized  polynomials  which  are 

t 0 

generated  by  a recursion  relation  identical  to  the  one  obeyed  by  the 
shell  states.  This  ortho-normalization  is  in  the  following  sense: 


(2.  16) 


Let  us  now  return  to  the  case  of  a composite  band  (such  as 
a 5 -fold  dt  band  complex),  the  single  particle  energy  £ and  the 
electron  charge  density  are  given,  respectively,  by 


<2.  17) 


0nir)=2,  t/u) 

D< 


(2.  18) 
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An  important  remark  about  this  method  should  be  made: 

The  self-consistent  Hamiltonian  must  have  the  property  that  the  total 
electronic  charge,  calculated  with  a unique  Ep  , precisely  cancel 
the  total  ionic  charge. 

With  this  self-consistent  £ and  'Yl(JC)  » we  can  calculate 
the  ground  state  energy  of  the  system  by  the  KS  equation: 


E=  £ — j Yp  ^X  + [ 'vujo  [ *1  <2. 19) 
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Chapter  III 


Example:  Stacking  Fault  Energy  of  a Nickel  Crystal 
in  3d^®4s®  Configuration 

In  this  chapter  we  apply  the  method  developed  in  Chapter  II 

10  0 

to  the  calculation  of  the  stacking  fault  energy  of  nickel  in  the  3d  4s 
configuration,  the  hypothetical  case  of  filled  <k  bands  (the  actual 
number  of  eJL  electron  in  metallic  Nc  is  about  9.4/atom).  Before 
proceeding  with  the  calculation,  let  us  briefly  review  the  properties 
of  a stacking  fault. 

A compact  structure  can  be  described  as  a stacking  of  hexa- 
gonal close  packed  planes  A,  B and  C;  f.  c.  c.  and  h.  c.  p.  lattices 
being  respectively  defined  by  periodic  successions  ABCABC  and 
ABAB  as  shown  in  Fig.  2. 

In  a f.  c.  c.  lattice,  there  are  two  types  of  stacking  faults. 

The  first  is  an  intrinsic  stacking  fault,  which  is  obtained  by  removing 
a plane  of  atoms,  i.  e.  ABCA1CAB  as  shown  in  Fig.  2(c).  The  second 
is  an  extrinsic  stacking  fault,  which  is  obtained  by  inserting  an  extra 
plane  C between  plane  A and  plane  B,  i.  e.  ABCACBCA  as  shown  in 
Fig.  2(d).  All  the  above  faults  preserve  close  packing,  so  that  any 
atom  keeps  the  same  number  of  the  nearest  neighbors  at  tlx  iame 
distances.  Thus,  it  is  reasonable  to  assume  that  the  formation  of 
these  faults  will  require  a small  energy.  From  the  experimentally 
determined  stacking  fault  energies  , they  are  seen  to  be  indeed 

small  compared  to  the  surface  energies  and  grain  boundary 


' 

: 


! 
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\z 


energies  ^ in  f.  c.  c.  lattices  (for  example,  in  Nl*.  Jjp  — 23  ° 

13  y 14- 

).  Although  they  are 

small  energies,  the  stacking  fault  energies  are  very  important  from 

a metallurgical  point  of  view,  since  they  govern  some  of  the  plastic 

properties  of  metals  and  are  useful  for  the  study  of  the  relative 

15 

stability  of  f.  c.  c.  and  h.  c.  p.  phases. 

Let  us  now  return  to  our  calculation.  We  want  to  obtain  a 
reasonably  accurate  expansion  for  the  atomic  wave  functions  'P,(£> 

( is  the  index  of  state,  i.  e.  s,  p,  d,  etc.  ) which  are  calculated 

by  Herman  and  Skillman^  in  terms  of  a number  of  Gaussian  type 
orbitals  : 


(3.  1) 


where  Ofl  denotes  the  types  of  the  Gaussian,  for  example, 
equals  1,  2 and  3 denote  s,  p and  d types,  respectively. 

H 


ii*l 


(3.  2) 


The  2M  parameters  and  Cj,  are  to  be  chosen  so  as  to  make 
the  approximation  as  good  as  possible.  The  most  convenient  criterion 
to  use  is  to  minin  ize  a weighted  squared  deviation, 

M 


D*Z.[u„trr)-Z. 


(3.3) 
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where  UOfp)is  a weight  function,  and  its  choice  is  governed  by  the 
use  to  which  the  expansion  is  to  be  put.  For  example,  our  purpose  of 
the  expansion  is  to  facilitate  a three-dimensional  integration  in  which 
Y is  the  distance  from  a given  center,  then,  in  Eq.  (3.  3),  we  would 
choo  se 

tocr)*rz  (3.4) 

which  is  proportional  to  the  volume  of  an  element  of  spherical  shell 
of  radius  Y . 

The  selection  of  the  best  set  of  the  exponential  parameters 
and  the  coefficients  presents  considerable  difficulties. 

This  is  due  to  the  existence  of  multiple  minima  for  the  value  of  D 
in  the  space  spanned  by  variational  parameters  { eli  > Cc]  . Therefore, 
it  is  not  claimed  that  we  have  obtained  and  listed  true  optimum  values 
of  [olCjCi]  in  table  1.  The  sets  in  table  1 are  actually  chosen  in 
such  a way,  that  they  fit  the  function  uniformly  in  the  whole 

range  of  T . 

The  minimization  is  performed  using  an  efficient  IBM 
Fortran  program  and  the  "best"  set  of  are  given  in  table  1. 

We  use  the  same  routine  to  expand  the  atomic  effective 
potential  in  terms  of  a number  of  "Gaussian"  which  are  shown  in 
Eq.  (2.  6)  and  the  exponential  parameters  and  the  coefficients  are 
given  in  table  2.  The  values  of  ^ for  a perfect  crystal  and 
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a stacking  fault  are  0.  148669087  au  and  0.  148669870  au,  respectively. 


In  the  case  of  a filled  band  the  integral  of  Eq.  (2.  14)  van- 
ishes by  Eq.  (2.  16),  except  when  ^=0  . Thus  the  expression  for 

'flu)  reduces  to  a trace  over  all  the  GWF : 


«2i2*  I ftjijnl2" 

i JL 


(3.  5) 


and  the  expression  for  6 reduces  to  a trace  of  Hamiltonian  over 
all  the  GWF: 


€ =2j2j  k l m i d,  y y (3  6) 

* 1 

where  <«*,  i I HU,  i>  is  the  self-energy  of  type  ck  at  site  £ . 

We  shown  in  Appendix  B that  the  stacking  fault  energy  for  a 

2 

filled  band  nickel  is  1.  98  ergs/cm  which  is  very  small  compared  to 

2 

the  measured  stacking  fault  energy  220  ergs /cm  . The  difference 
between  our  theoretical  estimate  and  the  measured  value  is  mainly 
due  to  the  fact  that  our  model  treats  Ni  ar6  a filled  band  metal.  We 
note  that  a vanishing  result  was  obtained  by  other  theoretical  calcu- 
lations^ using  a full  d-band. 


Chapter  IV 


9. 4 0.6 

Stacking  Fault  Energy  of  Nickel  in  the  3d  4s  Configuration 


In  this  chapter  we  apply  the  method  explained  in  Chapter  II 
to  a model  for  paramagnetic  nickel. 

We  know  from  the  energy  band  calculations^  that  there  are 
9.4  d electrons  per  atom  in  a solid  nickel  crystal.  Therefore,  we  will 
calculate  the  stacking  fault  energy  of  solid  nickel  in  the  3d^’^4s^'  ^ 
configuration.  In  this  calculation,  we  neglect  the  effects  of  sd  band 
hybridization  and  study  only  the  role  of  the  d electrons  in  the  stacking 
fault  energy  in  the  Hartree  approximation.  To  do  this  calculation  we 
follow  the  same  procedure  as  we  mentioned  in  Chapter  II  except  for 
constructing  the  atomic  effective  potential  V'  . Since  the  4s  wave 
function  is  highly  extended  in  space,  it  is  necessary  to  sum  several 
shells  ( > 10)  of  neighbors  in  order  to  obtain  an  accurate  value  of 
(see  Eq.  (2.  1)).  In  order  to  reduce  computation  time,  we 
separate  the  atomic  electronic  density  into  what  we  will  call  a local- 
ized atomic  charge  density  Tl^tt)  which  receives  contributions  from 
the  Is,  2s,  2p,  3s,  3p,  and  3d  states  and  a non-localized  atomic 
charge  density  T\.  (£)  which  is  made  up  from  the  4s  state: 


'n°‘in=  orctn  + ^ cn 


(4.  1) 


and  construct  V...  (£)  and  A V,ii  (.£)  in  the  following  fashion: 
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and  where^K^fJ^)  , yUj(Jl)  and  J^ll)  are  the  second,  third  and 
fourth  moments  of  the  site  1 , respectively.  In  this  scheme  w< 


we  use 


jc  g 

Eq.  (2.  13)  to  determine  the  Fermi  levels  Ep  and  Ep  of  the 
crystal  with  and  without  a stacking  fault.  The  stacking  fault  energy 


is  given  by 


r ( EF 

AE  = £\  E nf  (E)dE  -Z  E7lflE)<tE 

1 j L l J > 


= j Pce-e^[ 


_.  5F  - B 

2.0fl 

JL  JL  1 l 


M 


(4.5) 


where  the  superscripts  Sp  and  8 denote  the  system  with  and 

without  the  stacking  fault  energy,  respectively.  The  stacking  fault 

2 

energy  calculated  by  this  scheme  is  18.50  ergs/cm  . 

2.  The  non-self-consistent  scheme 

This  scheme  is  analogous  to  the  moment  scheme.  The  impor- 
tant feature  here  is  the  continued  fraction  representation;  it  is  one  of 
the  best  method  of  reconstructing  a density  of  states  from  a set  of 
continued  fraction  coefficients  (see  Fig.  3(a)).  This  reconstructed 
density  of  states  has  the  property  of  having  the  correct  leading  2 

moments  if  the  correct  leading  continued  fraction  coefficients  , 

bj  • are  used.  Using  the  density  of  state  'Ylg(E)which 

is  constructed  by  this  recursion  method  correct  to  ten  continued 

fraction  coefficients  and  Eq.  (4.  5),  we  obtain  for  the  stacking  fault 

2 

energy  a value  of  60.  20  ergs /cm  . 


-17- 


3.  The  self-consistent  scheme 


! 


We  know  that  in  a metal  the  perturbing  potential  due  to  a stack- 
ing fault  is  completely  screened  out  over  a few  interatomic  distances. 
The  Fermi  level  is  the  same  as  for  the  perfect  crystal.  However  in 
the  moment  and  the  non- self- consistent  schemes,  we  introduce  a non- 
physical variation  of  the  Fermi  level  for  a system  with  a stacking 

fault  to  satisfy  the  charge  neutrality  requirement.  In  order  to  avoid 

1 8 

this  problem,  we  use  a self-consistent  scheme  introduced  by  Allan 
to  insure  both  that  the  Fermi  level  of  a system  with  the  stacking  fault 
is  the  same  as  that  of  a perfect  crystal  and  that  charge  neutrality  is 
maintained.  Since  the  perturbing  potential  is  well  localized  near  the 
stacking  fault,  the  main  difference  between  the  self-consistent  and 
non- self-consistent  Hamiltonian  will  be  its  diagonal  elements 
< A I H \ <*>  x > where  X denotes  a stacking  fault  site.  We 

denote  this  difference  by  U.0  and  adjust  it  so  that,  using  the  bulk 

_B 

Fermi  level  fcp  throughout  the  stacking  fault  region  is  electrically 
neutral: 


**  IE,  lO  Xe  »£ 

1 l 


IE)AE 


(4.  6) 


In  the  presence  of  a stacking  fault  we  know  from  the  non-self-con- 
sistent calculation  that  should  be  a very  small  quantity  compared 

with  the  bandwidth.  Therefore,  we  can  assume  in  this  range  that  the 


j 


1 
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change  of  total  charge  at  the  stacking  fault  planes  is  linearly  propor- 
tional to  the  change  in  U.0  . Using  linear  interpolation  instead  of  the 
truly  self-consistent  loop,  we  obtain  UL0  = hi3>  *10  CL.U..  . 

The  stacking  fault  energy  in  the  self-consistent  Hartree 
approximation  is  given  by 

= Z \ E-Uo^c\e  - 2 \ FE/nf(E)o(E 

l J JL  ^ i J £ 

*-•5  | L (4.7) 

where  O'lo  and  \J0  are  the  electron  density  and  the  potential  in  the 
perfect  crystal.  In  the  crystal  with  a stacking  fault  these  quantities 
are  changed  into  and  V&+  <o\J  . The  first  two  of  Eq.  (4.7) 

corresponds  to  the  one -electron  contribution.  The  last  two  terms  of 
Eq.  (4.7)  are  added  to  prevent  the  double  counting  of  Coulomb  inter- 
actions occurring  in  the  integral.  It  can  be  written  as 

--5 1 j + ^Y)(£)*\/(£>&]  (4>8) 

It  is  easy  to  show  that  the  first  two  terms  in  Eq.  (4.  8)  give  an  equal 
contribution.  Here  we  must  notice  that  is  well  localized  near  the 
stacking  fault  and  the  preponderant  term  in  'Vl0  is  ^=0  (see  Eq.  (2.  14)). 
Therefore  we  can  rewrite  Eq.  (4.  8)  under  the  form 

— ?b~)U.o  (4.9) 

where  and  ?<sn  are  the  number  of  d-electrons  per  atom  in  the 
bulk  and  on  the  stacking  fault,  respectively.  The  stacking  fault 
energy  calculation  by  this  scheme  is  64.  28  ergs/cm  . 


Chapter  V 

Comparison  and  Discussion 

For  purposes  of  comparison  we  have  drawn  the  bulk  density  of 
states  computed  with  ten  continued  fraction  coefficients  (a.  b-- 

j =0,  ....  9)  and  a histogram  deduced  from  Pettifor19  with  the  appro- 
priate scaling  in  Fig.  3(a).  The  good  agreement  of  our  results  for  the 
bulk  density  of  states  with  band  structure  calculations  gives  us  confi- 
dence in  our  further  discussion. 

We  have  also  calculated  the  stacking  fault  energies  with  these 
three  schemes  which  were  discussed  in  Chapter  IV  using  the  Slater- 
Koster  overlap  parameters  <U<r  , U*  and  «W$  to  calculate  the 
nearest  neighbor  hopping  integrals  instead  of  using  our  constructed 
Hamiltonian.  The  parameters  we  have  used  are  (in  atomic  units) 

(US's  -0.020238  <Utt  = o,o0<}42>  d4$  = -o,  oollf>l$ 

which  was  chosen  by  Desjonqueres  et  al.  20  The  results  of  these  cal- 
culations are  given  in  Table  3. 

From  Table  3,  we  can  draw  the  following  conclusions: 

1.  To  replace  our  constructed  Hamiltonian  by  a parameterized 
Hamiltonian  seems  to  have  very  little  effect  on  the  stacking  fault 
energy. 

2.  The  density  of  states  which  is  correct  to  the  fourth  moment 
is  not  good  enough  for  the  stacking  fault  energy  calculations. 
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3.  The  difference  between  the  numerical  values  of  the  stack- 
ing fault  energies  obtained  in  the  non- self-consistent  scheme  and  that 
obtained  in  the  self-consistent  scheme  is  very  small. 

4.  The  results  we  find  are  smaller  than  the  available  experi- 

2 2 

mental  values  (the  range  of  values  = 1 6 0 ergs/cm  ~ 300  ergs/cm  ). 
Two  possible  sources  of  the  discrepancy  are  the  following:  our  neg- 
lect of  s-d  hybridization  and  the  fact  that  the  measured  values  refer 
to  nickel  in  its  ferromagnetic  state  rather  than  the  paramagnetic  state 
which  we  used  in  our  calculation. 

5.  Although  our  value  for  the  stacking  fault  energy  is  rather 
smaller  than  the  best  experimental  value,  it  is  in  much  closer  agree- 
ment than  is  the  value  calculated  by  Cyrot- Lackmann  using  the 
moment  scheme.  Our  result  is  a factor  of  three  higher  than  her 
result. 

Therefore,  the  next  step  in  such  calculations  would  be  to 
include  s-d  hybrdization  effect  for  a ferromagnetic  state  nickel. 

Let  us  finally  notice  that  our  scheme  could  in  a straight- 
forward way  be  used  to  calculate  energies  of  more  complicated 
defects,  for  instance,  formation  energy  of  vacancies,  etc. 


, 


\ 
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Appendix  A 


r 


In  this  appendix  we  present  expressions  for  the  integrals  used 
in  the  computation  of  the  matrix  elements  of  the  Hamiltonian. 

In  the  Cartesian  coordinates  a Gaussian  type  orbital  (GTO) 
centered  at  A is  written  as 


. (a. i) 


are 


Here  ^ A and  are  the  components  of  a position  vector 

ffl  relative  to  A , Ts  = T - A , and  Jt  , on  and  Ofl 

c*  ^ r*, 

"quantum  numbers"  and  are  zeroes  or  positive  integers. 

First,  let  us  define  some  notations: 

^ U _ i ^ ^ • ">  « ~ ««  N 


= ^+M*r'6cf+pBxr  = ^ (j,  u,  p^pe^f 

a 


(A.  2) 


where  PAX  and  PBX  are  the  component  of  the  vectors 
( PA  = JP  ~ A)  and  PB  [ (^B=  £ -B)  , respectively. 

I 


?A 


u^JPX|>  l-tal)Au.  . 


(A.  3} 


Then  we  have  the  following  expressions: 
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S ( 'nO— J A,,  (. S/,,  'McU^ 


34  - AB1 


X PA)J  7 PV  U«** 


^l'«ltM2')^ 

*ZIZ 

■^=0 


(A.  4) 


where 


^ = o<|  + cV2 

P = w>A+rf,g)/j' 


(*i-0!!  = l>3-5 lsi-0 


(A.  5) 


and  where  w mean8  "largest  integer  less  than  or  equal  to  OC  . 


r 


^' )] /•tM'AvwtC 


- 2 0<2  ^ S ( l*  + 2;  ^2^z)  + S ( ij/,/Wa+  2,  + S (.  l»;  <*,,  012+  2)] 

- 1 \ t(  t-l)  $ ( k-Z,™!,  ^z)  + 'Wal^a-O  S(  la,'1"*-*,  ^ 

= ?e  * 

iM 

*2  A>a,v(0%/*%  V 


(A.  6) 


cp* 


(A.  7) 


it, to 
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where  J)  — A.  'j  -f  'fe  - 2 (T+A+  A.}  — ^ U.+  'V'-t  to)  ^ 


and  where 


Ai^r,  u.  ^*3 


= C-oHt  (!,,!,, PA^PS,) 


r-  - . (-0 *i!  PCx 


L-2T-2U, 


( At) 


r*u_ 


r\  a i U-^r-2aV. 


(A.  8) 


In  Eq.  (A.  7 ) the  summations  with  respect  to  indices  A , ^ and 

extend  from  0 to  H,+  l,  . U/s)  and  £ U-20/al  , respec- 

tively. The  range  of  Ij.A/U 0 or  i*,t;  co)  can  be  easily  found  in 
the  same  way. 


fjr  jr 

e 3b,  Ob,  |r2_r, 

rv  i~*  ' 


*Z  B, 


S*-«Ar«ya 


u ^0  ^2;  Us,®4,^<,^xW 
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* Z B*  4 u ^ W» ' W*  WO&  ($ ) 


(A.  9) 


where  Jj  s **|+*at  MAi+MO-U-TT-IO 


£|  -«'|t  o(i 

p _ 

~ " fr, 

J-  JL+  J_ 
6 4r,  4^ 


£a  = °^5  -t  °<4- 

Q _ ^ieHE 
~ ^ 


(A.  10) 


and  where 


^l,,L„ r„r„t/  *'j  h,  Ax,Bx,Px,Sf,|  l5,U/x,px,a^f:) 

= (-if1  fi((  »v  Mx  ,PBx)  \i}  U!;U>  QCx , «P*) 

x i-.!  lx!  M,j‘  Kfrfi  Xltlw 

WO1'  (W*^1*  V.r,!  (i,-2r,)!  U»-*W! 


r ti4Aa-2(r,+r«.,)-2a  cu 

xfWviraiM!  =f 

J u.»  \ l,4  jL2-2ir,+r»')-^u]  i 


(A.  11) 


-26- 


I 


In  Eq.  (A.  9)  the  summations  with  respect  to  the  indices  i,  , ^ 

ra  , and  U.  extend  from  0 to  A,+  A,  , JL5+  ^ , [ L./a-]  , £ U/2~l 
and  C Ui+  iz)/z  - !Tt -V*"]  , respectively.  The  range  of 

•*»  ■*»/*)  or  to)  can  be  found  in  the  same  way. 


Appendix  B 


The  stacking  fault  energy  of  a nickel  crystal  is  the  difference 
between  the  ground  state  energy  of  a system  with  a stacking  fault  and 
that  of  a perfect  system: 

= ESF  -ES 


e*  - e6 


- A f ^ tP  ir  i,'  + A f *Q*V)  Jr  J,' 

"J  Ift'l  ~~  J ll-f' 

+ j ^tn[  ‘h 


where 


= ^ E,  -V  ^ Ea  + ^ Ei 


^E,s  t - 6 


(B.  1) 


^E5-  j ) [ 6 « , 1^11 )) -y#* ill)] k -J ,n.V)[f kI,11' 


(B.  2) 
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Here  the  subscripts  SF  and  B denote  the  system  with  the  stacking 

fault  and  the  perfect  system  respectively. 

10  0 . 

For  a nickel  crystal  in  the  3d  4s  configuration  the  single 

I] 

particle  energy  6 is  given  by  Eq.  (3.6).  Then  ^E,  is  given  by 

1.1?  M¥/v»? 

Since  ^ Ea  involve s two  six-dimensional  integrals,  usually 
it  is  impossible  to  do  this  integration  numerically.  However,  in  our 
representation  the  charge  density  OfUf)  is  expressed  in  terms  of  a 
number  of  GTO  (see  Eq.  (3.  5))  and  this  integral  can  therefore  be  cal- 
culated  analytically.  Due  to  the  large  computation  time  (because  the 
number  of  GTO  is  too  large)  we  cannot  in  practice  calculate  ^E2 
directly.  Here  we  want  to  give  some  estimates  for  this  term  under 
the  assumption  that  the  GWF  do  not  overlap  between  the  third  nearest 
neighbors.  If  we  write  (SB  , the  atomic  wave  functions 

can  be  written  as 

yn =<(£-!) 

The  GWF  and  the  electron  charge  density  are  given  by 

yt)=yu-^yt>(%'.y  , 


(B.  3) 


(B.  4) 
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Therefore,  the  coulomb  interaction  energy  can  be  written  as 


f j , 2 , t„u.'T»w  't'~  " " jrj.< 

J lt-1'l  ' A di- "t 

^RittlteWte  .JrJ/ 


*2.2.2.  Wf/ft) 

(Si*'  7 

S*Jr' 


Ir-r'l 


(B.  6! 


The  first  term  on  the  right  hand  side  of  Eq.  (B.  6)  is  a two-center 
integral  ( ($,!{)  and  the  second  term  is  a three -center  integral 
( ft  )•  Since  the  stacking  fault  preserves  the  same  number  of 

nearest  and  second  nearest  neighbor’s  at  the  same  distance  as  that  in 
perfect  lattice,  the  two-center  integrals  give  the  same  values  for  the 
system  with  and  without  a stacking  fault  (see  Eq.  (A.  9)).  By  the 
same  reason,  the  three-center  integrals  give  the  same  results  for 
these  two  different  systems  except  in  the  following  situation:  the 
three  centers  involve  three  parallel  planes  in  which  the  stacking  fault 

x 

occurs.  The  difference  of  this  kind  of  term  between  the  system  with 
a stacking  fault  and  the  system  without  a stacking  fault  is  bounded  by 
0.  16  ergs/cm  . Therefore,  we  can  neglect  the  effect  of  on  the 

stacking  fault  energy. 
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Figure  1.  Flow  chart  for  the  calculation  of  Ni  stacking  fault  energy. 

Starting  from  atomic  wave  function  and  densities,  we  first 

»«( 

obtain  (a)  generalized  Wannier  function  U.^(.r)  and  (b)  an 
effective  potential  for  the  stacking  fault.  From  these  the 
Hamiltonian  matrix  and  finally  single  particle  energy,  elec- 
tron density  and  total  energy  are  calculated. 
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Figure  2(a).  Projection  normal  to  the  (1,  1,  1)  plane  showing  the 
three  types  of  stacking  positions  A,  B,  and  C;  (b)  Perfect 
f,  c.  c.  stacking;  (c)  Intrinsic  stacking  fault  in  f.  c.  c.  struc- 
ture; (d)  Extrinsic  stacking  fault  in  f.  c.  c.  structure. 


-36- 


Figure  3 (a).  The  bulk  density  of  states  of  Ni.  The  histogram 
is  deduced  from  Pettifors  exact  calculations,  the  full  line 
being  our  result  correct  to  the  lol^  continued  fraction 
coefficients. 
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Figure  3 (b).  The  local  density  of  states  on  the  stacking  fault 
with  U0=O 
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Figure  4.  The  bulk  density  of  states  given  by  the  Edgeworth  series 
fitted  with  4 moments  for  Ni. 


Table  1.  Orbital  exponential  parameters  and  expansion  coeffi- 

cients Ci  for  Is,  2s,  3s,  2p,  3p  and  3d  orbitals  in 
atomic  unit 


oLi 

Ci 

8 

1 s 

2s 

3 s 

11516.  5 

33. 0481 

8. 82768 

4. 20507 

3061. 22 

53. 2612 

19.3895 

5. 94182 

962. 010 

40. 2191 

6.74675 

6.  48046 

530.  057 

64. 2795 

32.  1637 

7. 55527 

165. 397 

80.  7803 

4. 20110 

19. 2241 

134. 438 

-19. 0622 

28.  9506 

-6. 32380 

59.  5449 

14.  3051 

-.969386 

3. 35939 

15.  1521 

-.  259403 

-14.  7299 

-12.  5191 

6.  24163 

. 165880 

-3.  28423 

-. 048757 

1.  39054 

-. 033213 

-.  028432 

3. 56488 

P 

2p 

3p 

340.  262 

116.  178 

46.  2843 

155.  768 

154.  761 

59.  3034 

50.  4147 

127.  228 

55. 0839 

20.7513 

62.  4323 

26.  2800 

7.  06297 

9.45597 

-.  526122 

2.  92262 

-. 394319 

-7. 10654 

.951675 

. 035764 

-1.  23424 

d 

3d 

48.  9403 

62.747 

13.7169 

36.  546 

4.  63951 

13. 677 

1.  57433 

2.  6311 

.486409 

. 24164 
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Table  2.  Parameters  of  atomic  effective  potential  expansion  in  atomic 
unit.  o(i  are  the  exponential  parameters;  Ci  are  the 
coefficients 


( 

1. 710018091851 
26.72930960088 
7. 640172520328 
2. 324012025830 
0. 260972626075 


Cc 

-26. 39857137242 
17.  08366865139 
21.  12143227226 
32.  21210656087 
-0. 861105554065 


z 

Table  3.  Comparison  between  the  stacking  fault  energies  in  ergs /cm 
which  are  obtained  by  the  moments  scheme,  the  non-self- 
consistent  scheme  and  the  self-consistent  scheme  by  using 
(a)  our  constructed  non-periodic  Hamiltonian,  (b)  3 parame- 
ters d.A.0-  , dLAtr  and  AA5  fitted  Hamiltonians. 


Moments  N.  S.  C.  S.  C. 


a 


18.  50 


60.  20 


64.  27 


